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In this talk, I will discuss several issues related to the theory of the elec- 
troweak phase transition. A discussion of other aspects which we have recently 
addressed^ '^ can be found in the article by A. Linde in these proceedings. 

The theory of phase transitions in the early Universe was pioneered in the 
seminal work of Kirzhnits and Linde in 1972. Early works claimed that the phase 
transition corresponding to the breakdown of SU(2) x U(l) was of the second 
order. '' Finally in 1976, it was shown to be of the first order. This fact essen- 
tially lay dormant for many years until it was realized that this character of the 
phase transition can provide the non-equilibrium environment essential for produc- 
ing a baryon asymmetry at these scales; the baryon violating process is typified (at 
temperatures below the critical temperature) by the sphaleron. We now believe 
that the phase transition proceeds by the nucleation of critical true vacuum bub- 
bles from a rather weakly supercooled state. These bubbles expand under pressure 
forces and quickly fill space. Baryogenesis typically proceeds in the bubble walls 
where is changing,^" ^^ or in a thin layer in front of the wall.^^ 

It is potentially important for a detailed understanding of the baryon asym- 
metry, to understand the propagation of the bubble wall, namely its size and shape, 
and its terminal velocity (if any). In the last half of this presentation, I will discuss 
our estimates of these quantities, and compare to other recent analyses. Work is 
continuing on refining the analysis and considering new phenomena that may be 
important. ^^ 

In the first section of this presentation, I will review the standard analy- 
sis leading to the derivation of the effective potential for the Higgs mode. It has 
recently been emphasized that higher order corrections are potentially im- 

portant because of infrared divergences and indeed can lead to the breakdown of 
perturbation theory in relevant regions of parameter space. We discuss several ap- 
proaches to resuming perturbation theory that soften the poor infrared behaviour. 
Finally, we outline by powercounting the region of parameter space for which the 
resummed perturbation theory is good, and comment on what effects the analysis 
has on various models of baryogenesis at the electroweak scale. 

1. The Electroweak Potential 

1. 1. One-Loop Analysis 

The standard approach to equilibrium properties of finite temperature field 
theories is the imaginary time formalism. In this technique, one continues to 
Euclidean spacetime, and compactifies the time coordinate on a circle of radius 
U(3 = T~^. The resulting path integral is then identical to a Boltzmann sum over 
states. The fields in the theory may be expanded in Fourier series over a discrete 



set of Matsubara frequencies 

^(f,r) = J]e^"">n(f). (1.1) 

n 

Requiring the fields to satisfy (anti-)periodic boundary conditions on the circle 
leads to uJn = (2n + l)7rT (fermions) or uj^ = 27rnT (bosons). In this way, the 
theory can be thought of as an infinite set of coupled three-dimensional fields with 
a hierarchy of masses m^ = u;^ + m^ where mo is the usual zero temperature mass. 
As we will see later, this viewpoint will be of particular value to us. The quadratic 
part of the action (for bosons) may be written as 

Sl3 = iTj2 I d?x (t)^n{-V^ + ool + mt) (t)n. (1.2) 

The effective potential for the Higgs field can be computed by putting 
in a source and studying the response of the theory to this source. At the one 
loop level, this is equivalent to shifting by its saddle point value; the familiar 
result is just the determinant of quadratic fiuctuations. That is, at lowest order, 
we can neglect interactions and the effective potential is that of a collection of free 
particles interacting with a heat bath at temperature T with masses given by 0. 
We have 



Vr,i-/oop(0) = Vb(0) - iT^lndet [-V^ + o;^ + ^2(0) 
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where 



oo 

I^{y) = ± f dxx'^lnflT e-^y^^'+y'] . (1.4) 


Expanding for small y, i.e., m{(f)) <^ T, we find the following potential: 

\/t(0) = D{T^ - T^o)^'' - ET<i? + iAr04 (1.5) 

where D, E and \t are determined in terms of the gauge and Yukawa couplings 
of the theory and the Higgs self-coupling. We will assume throughout that the 
Higgs boson makes only a very small contribution to the potential, i.e., that it is 



sufficiently light. The parameters in the effective potential, Eq. (1.5), are given 
in many previous works, including Refs. 1,2. For present purposes, we note the 
following. The parameter E is crucial in making the transition of the ffist order; 
when E is non-zero there is a second minimum at positive which first appears at 
a temperature Ti, becomes stable at a temperature Tg and then becomes the true 
vacuum below the temperature Tq. Thus, we expect that supercooling can take 
place as the Universe cools, the field getting stuck in the false vacuum state at 
= for T < Tc. In a second order transition, the parameter E is zero, and there 
is a smooth transition from = to 7^ with no supercooling. In our context the 
supercooling is crucial for supplying the non-equilibrium environment necessary for 
baryogenesis. The transition from false vacuum to true vacuum proceeds by the 
nucleation of bubbles of true vacuum which then expand (if large enough) to fill 
the universe, baryons being produced in the vicinity of the advancing bubble walls. 
In the minimal standard model, the parameter E, taking into account the 
analysis presented above, is given by 

i? = ^(2m3, + m|). (1.6) 

Only bosonic fields contribute. This is understood in the present context by the 
appearance of a non-analyticity in the dependence of /_|_ on y'^ coming from the 
minus sign preceding the exponential in Eq. (1.4). This non-analyticity is absent 
for J_. It will be of some importance to us to understand the origins of this non- 
analyticity. As I will now demonstrate, it is due to infrared divergences in loop 
integrals. 

The cubic term in the potential may be derived in several other ways by 
evaluating Feynman graphs. In order to do this, we will find it convenient to 
evaluate, instead of the potential directly, the tadpole graphs that contribute to 
dVT/d(f). In the minimal standard model at one loop, we then find 
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where the hi are the couplings of to the various particles and Pi is the Euclidean 
fermion propagator. The trace gives a factor of 'mi{(f)) = hj^icj) for the fermion 



contribution. Re-writing hi = hr^ for the fermions, we obtain 
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I have separated off the a; = terms as these have different behaviour for small 0. 
Indeed we can perform the above integrations by dimensional regularization and 
the frequency sums can be zeta-function regulated to obtain, for small mi{(f))/7iT, 
the following: 

^ = -^ E hiT<j)mi{<P) + J2cih^T^<l^+■■■ (1-9) 

i, bosons i 

where the q are 0-independent constants. We note that the first term, which came 
entirely from the zero frequency mode of the bosons, integrates to give a cubic 
term in Vy, whereas all of the other modes give corrections to quadratic (and 
higher) terms in 0, and in particular, are analytic in 0^/T^. Thus the cubic term 
is understood as arising from the infrared region of 4-momentum space. 

Now, if we go on to powercount higher loop graphs, we find that the infrared 
behaviour found at the one loop level worsens. Indeed, one finds that there is an 
effective expansion parameter that is of order g T /m {(p) ~ (T/0) , where g is 
some generic coupling constant; for concreteness, we take it to be the SU(2) gauge 
coupling. The appearance of this expansion parameter means that perturbation 
theory breaks down at less than or of the order of T. This is clearly a disaster 
for an examination of the phase transition, where everything typically happens 
in this regime. ^"^ In what follows, we will consider improvements of perturbation 
theory obtained by attempting to take into account all of these infrared divergent 
contributions in some consistent way. To begin, we will assume that the coupling 
g is sufficiently small that perturbation theory using it as an expansion parameter 
makes sense. In contrast to ordinary 4-dimensional field theory, we will see that 
one does not here have the luxury of factors of 47r in the expansion parameter. 
For these reasons, much of the powercounting arguments supplied below will be 
merely formal arguments valid for small g\ an accurate numerical analysis awaits 
the future. 



-k In particular, it has been claimed that (light) Higgs bosons may induce numerically large 
corrections. See Refs. 20. 



1. 2. Infrared Improvements 

The fact that higher loop graphs are important near the phase transition 
was first noted by Weinberg"^ and was also discussed in Refs. 6,22. Indeed in 
order that a theory develop a second minimum at low temperatures, it must be 
that the perturbative expansion breaks down because the second minimum is 'non- 
perturbatively far away' from the original vacuum state. More recently, this has 
been emphasized in the context of the electroweak theory by several authors. In 
Ref. 17, Brahm and Hsu attempted to take into account some of the higher order 
graphs but found a large negative linear term ~ —g^T^cj), leading to the conclu- 
sion that the transition is at most very weakly first order. In a similar analysis, 
Shaposhnikov found a large positive linear term ~ g T (p, leading him to con- 
clude that the transition could be very strongly first order. In work coincident with 
our own, Carrington^^ found corrections to the cubic term, leading her to initially 
conclude that the transition is more strongly first order. All of these results have 
now been understood to be incorrect, as I will explain in detail below. In particular 
I will show that linear terms never arise, but that the cubic term is modified in an 
important way. 

Let us now look at higher loop graphs in more detail. The graphs of Fig. 1 
where the large loop is at cj = and the smaller loops are at a; 7^ are the most 
infrared divergent graphs. These have been referred to in the literature as 'daisy' 
or 'ring' diagrams. If we powercount the graphs for m{(f)) <^ ttT, we find that they 
are proportional to 

22^2 \P /2^\2p 
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^ring,p~ foTTT ~ T (1-10) 



and thus are not formally suppressed by powers of coupling constant. To proceed, 
we must resum perturbation theory in some way. In the following sections, I will 
discuss in detail two methods of understanding this resummation. 

1.2.1 Schwinger-Dyson Approach 

One notes that the diagrams in Fig. 1 correspond to insertions of one-loop 
propagator corrections on the large loop. This leads us to believe that perturbation 
theory may be improved if we can take into account these loop corrections to 
the propagators. In an idea outlined in the book by Kapusta,^^ and stressed by 
Carrington^^ in the context of the Standard Model, one extremizes the effective 
action, not just with respect to the vacuum value of the Higgs field, but also with 
respect to the polarization tensors, {Il{uJn,k)). Thus we go one step beyond the 



mean field approximation. To do this one writes a shifted propagator 



'^o,lik)= P + ujl + ml 
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where we for simphcity give exphcit formulae for a scalar boson only. The result 
can be immediately generalized, and indeed we will do so later for the Standard 
Model. We now add zero to the Lagrangian in the following way: 



C{VQ)=C{V)-\TY,(t>~n^r. 
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We are to do perturbation theory with the first term, and treat the last term as 
a counterterm. We can then compute the one-loop effective potential, which looks 
like2i 

r (fk °° 

(1.13) 

Note the important correction term to the one-loop logarithm. The integer i counts 
loops, and one should note that these higher loop diagrams involve corrected prop- 
agators. The self-energy is obtained self-consistently by requiring that it extremize 
the potential. An arbitrary variation of this effective potential can be written 
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where and 11 have been varied independently. Thus the requirement on the 
self-energy turns out to be just the Schwinger-Dyson equation 
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given diagrammatically in Fig. 2. 



The variation with respect to gives us the equation for the tadpole, shown 
in Fig. 3, where the internal lines are full propagators, and the vertices are uncor- 
rected. To show that the perturbative series has indeed been improved, we can 
evaluate, again in just the simple scalar theory (we return to the realistic case 
presently) this expression at one loop; we find 
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where we have isolated the zero-frequency contribution as above. Taking the limit 
A; ^ within 11, we obtain 
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and hence a 'cubic' term 
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As long as no(0) is not singular in the limit — ^ (it isn't), we obtain no linear 
terms in 0. Furthermore, if no(0) is large compared with m^ (true for some region 
of small 0), the cubic term is reduced to essentially zero. This will have important 
consequences later for the electroweak potential. In the scalar theory, Eq. (1.15) 
gives us at lowest order 



no(g ^ 0) = 3AT2 I J2 I 7^^-(^) 1 + 6A02p_„(g - k) 
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The terms in the ellipsis involve corrections of order m/T, and include dependence 
on Iln. The leading term, of order AT^, is unambiguous, and we see that the 
would-be cubic term is removed for small 0. 

We now consider the Standard Model. In this case, we run into additional 
problems: in a gauge theory, the 'magnetic mass' vanishes. That is, if we compute 
Iio{k -^ 0) for the transverse gauge bosons, we find there is no term of order g'^T'^. 
We will see in the following sections that this fact leads to an as yet insurmountable 
difficulty with the perturbative expansion, at least in some range of 0, and has a 
direct analogue in finite temperature QCD. 



Before discussing this further let us note the changes to the electroweak 
potential coming from the resummation. The top quark self-energy is Iltop — 
ggT /6, but being a fermion it has no cj = mode, and so does not contribute to 
the infrared problem. Higher order graphs lead to numerically small corrections to 
the 0^ and 0^ terms. To discuss the gauge bosons, let us work in Coulomb gauge. 
In this gauge, we have propagators that mix the transverse gauge fields with the 
Coulomb mode and the Goldstone mode. At a; = 0, these decouple and we find 
that the Coulomb and Goldstone mode self-energies have a term of order g'^T'^, as 
does the physical Higgs, but the two transverse modes have no such term. This is 
as we discussed above. Referring to the result of Eq. (1.17), we write the corrected 
cubic term from the W^ as 



K=ubic - ^^^T 



{ml. + mlf/^ + 2{ml.f'] (1.19) 

where rriD = y^^oik ^ 0) ~ gT is the Debye mass. There is also a similar 
contribution from the Z . As we will see later, this is reliable for ^ gT, and 
thus we see that the cubic term has been reduced by a factor of 2/3, as the first 
term in the brackets of Eq. (1.19) gives (small) corrections to the and terms. 
This has important consequences for baryogenesis. We note that the position of 
the minimum is now given by 



T 
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where e = {T^ — T"^) / {T^ — T^) gives the temperature of nucleation. The position of 
the minimum, for a given value of T has been moved in towards the origin and hence 
the transition is less strongly first order (the height of the barrier is correspondingly 
lower). The bound on the Higgs boson mass, found by requiring that sphalerons not 
wash out^^ the baryon asymmetry after the transition, is expected to be lowered 
by a similar factor, down to about 40 GeV. We conclude from this analysis that 
the minimal standard model is incapable of producing sufficient ub/s (regardless 
of the amount of CP violation) because of the LEP bound of ttih ^ 57 GeV. 

This bound is easily avoided however by extending the model. A multi- 
Higgs model (two or more doublets) is particularly attractive here, because it can 
easily avoid the mass bound^^'^^ and there is CP violation at a (nearly) sufficient 
level. Anderson and Hall^^ have also noted that adding singlet scalars can modify 
the dependence of the Higgs mass on the self-coupling in such a way that the 
bound is avoided. As well as raising the theoretical bound on the Higgs mass, an 
extended model does not suffer from such a restrictive laboratory bound as does 
the minimal standard model (because of, e.g., tan/? dependence, etc.). 



1.2.2 Infrared Effective Action 

It is also useful to understand these aspects of the electroweak potential in 
an alternative way. The idea here builds on the fact that the interesting behaviour 
comes from the infrared region. Since, at weak coupling, we are interested in energy 
scales much less than ttT, one can first integrate out all modes of the various fields 
with non-zero Matsubara frequency, thereby obtaining a three-dimensional effective 
action for the zero-frequency bosonic fields. As we will see below this technique 
automatically sums the troublesome graphs, and automatically gives the correct 
combinatorics; we are able to clearly see the absence of linear terms in the potential. 

Integrating out all of the heavy modes corresponds to evaluating all Feyn- 
man diagrams containing loops of heavy fields, with an arbitrary number of ex- 
ternal light legs. Thus we obtain a three dimensional effective action involving an 
infinite number of operators. In particular, the quadratic terms pick up important 
contributions. In a scalar theory, this action is given by 



SCq^ad = h j d^x(l)l 3AT 5^ y 
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The quadratic part of the action thus has a term of order AT and all corrections 
to this are analytic in (m/T)'^, because the n = contribution is absent. The 
contribution of the n = mode corresponds to loop graphs in the infrared effective 
theory. 

Having integrated out the heavy modes (at one loop level), we can now 
compute the one-loop effective potential in the effective theory. This is just the 
determinant of Gaussian fluctuations in the low energy theory, and corresponds to 
the bubble graph of Fig. 4; if one writes this out in terms of the heavy lines that 
have been integrated out (Fig. 4), we see that we have reproduced the sum of the 
dangerous ring diagrams. One can easily check that the cubic term is correctly 
reproduced, as in Eqs. (1.17) and (1.19). 

In the case of the Standard Model at one loop, one finds results similar 
to Eq. (1.20) for the Coulomb and scalar lines and thus one has the Debye mass 



* We have used dimensional regularization. 
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corrections of order g T for these fields, whereas the quadratic term for the trans- 
verse gauge bosons goes to zero with k and 0, as discussed above. These "mass 
corrections" are analytic in |0| , in particular, they do not contain linear terms in 

0, and thus there are no linear terms in the effective potential, at least at this order. 
In principal, we should go on to examine higher loop graphs. In the next section, 
we will discuss the improvements that have been made to perturbation theory by 
summing the ring diagrams. In the present context, however, we can easily under- 
stand the absence of the linear terms at the two-loop level as well. The dangerous 
graphs are those with a transverse gauge boson in one loop and some other field 
with a Debye mass in the other loop. Individual graphs are badly behaved, that 
is, they apparently give linear terms, but the sum of all such two-loop graphs is not 
so singular. This can be understood by thinking of the loop containing fields with 
a Debye mass as a one-loop correction to the propagator of the transverse gauge 
boson, i.e., by an insertion of the polarization tensor. This has been calculated 
at one-loop level and contains a momentum factor in the numerator which softens 
the infrared divergence of the two-loop graphs. Thus one sees that there are no 
linear terms in the effective potential through order g (pT . 

Additional higher loop graphs are potentially troublesome; these are known 
in the literature by the name "superdaisy" graphs. In the present context these 
are taken care of by further integrating out modes. Recall that we have integrated 
down to a momentum scale of order ttT. In this theory, we have fields with masses 
of order gT, and fields which are lighter. We can consider moving the cutoff lower, 
below gT, whereupon the Coulomb and Goldstone modes are integrated out. As 
long as the temperature is not too close to To, we can also integrate out the Higgs 
boson as well, leaving only the transverse modes. With such a low cutoff, the 
problems arising from the superdaisy graphs are moot. In the next section, we 
will study the infrared behaviour of the full theory and the validity of perturbation 
theory by looking at this low energy theory. At temperatures in the vicinity of the 
phase transition, additional divergences arise because the Higgs boson is becoming 
light. In what follows we assume that the Higgs boson mass is at least of order gT. 

1. 3. How good is perturbation theory now? 

We would now like to discuss the improvement of the perturbation expansion 
in a gauge theory and discuss whether or not we can reliably predict a first-order 
transition. We have noted that the expansion fails because of the lowest-order 
vanishing of the magnetic mass of the transverse gauge bosons. We expect that 
the associated infrared divergences will be cut off in some way. For the purpose 



f Recall that the loops in the effective theory involve only lo = Q modes; the loop momenta 
are cut off at the Debye mass. 
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of the present discussion, let us suppose that this cutoff is provided by a magnetic 
mass TTiniag ~ g'^T. Consider then the theory where we have integrated out all 
modes except for the transverse gauge bosons (i.e., at scales below gT). If we 
powercount vacuum graphs, we find that they behave like 



'1.21^ 



If we neglect rnmag, the expansion parameter is of order gT / (p, and the expansion 
is valid for ^ gT. This range of is consistent with neglecting my/ in relation 
to TTimag- We can interpret this in another way. In the vicinity of the phase 
transition, 0/T ~ E/\ ~ 5'"^/A, and so the expansion parameter is of order \/ g^. 
We conclude from this that large mjj is a potential problem for us, at least if we 
want the minimum of the potential to appear in the perturbative regime. 

The typical shape of the potential is shown in Fig. 5. The exact placement 
of the limit of perturbative applicability and the position of the minimum depend 
crucially on numerical factors. In any case, for sufficiently strong first order tran- 
sitions (where 0/T ^ 1) we can reliably predict the presence of a minimum away 
from the origin, and we expect a first order transition to occur. Further evidence 
is provided by the e-expansion and the absence of an infrared fixed point. ^^ 

We can also worry about the presence of bad behaviour at small 0, and 
whether or not symmetry is truly restored at high temperatures. Here, the be- 
haviour of higher order graphs is similar to that of QCD. There, the potential 
can be written, as in Eq. (1.21), as 



V = ... + T''{ 0{g'') + Oig" \ng') + O 



fY.^g'T/m^,,r 



where all of the last terms are formally of the same order. The standard wisdom is 
that we can compute the small 0{g'^) and 0{g^\ng'^) corrections, and then hope 
that all the rest adds up to something small. This is an unsolved problem. In the 
case of the electroweak theory, we arrive at the same result, but we have additional 
information. The cutoff at 0{g'^T) implies polynomial behaviour in for small 



I For temperatures very close to the transition, the Higgs wiU also become light, and the 
perturbation expansion will break down. We suppose here that we are at sufBciently high 
T that the Higgs may be integrated out as well. 
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and thus we expect, for example, the quadratic term can be written 



K 



quad 



-W + '^Bvl + DT^ + 0{gHng^)T^ + yq\g^T/mra.^Y 






Even though we cannot hope to compute past 0{g'^\\ig'^) here, we do not expect 
any bad behaviour at small 0, and thus we can claim that symmetry is restored 
at high temperature, regardless of infrared divergences, provided they are cutoff in 
some way. Furthermore, as above, we conclude that the transition is first order, 
as long as it is sufficiently strongly so. This includes those theories capable of 
producing a baryon asymmetry. 

2. Bubble Wall Propagation 

A detailed study of baryogenesis at the electroweak phase transition neces- 
sarily involves an understanding of how and when the bubbles are formed in the 
first order transition, as well as how they propagate. In this talk we will give only 
a cursory discussion of the former, noting only those aspects that are of relevance 
to the propagation of the bubble walls. 

Studies of the nucleation event^'^^ indicate that bubbles are formed at a 
temperature just below T^, but sufficiently far below this that the thin-wall ap- 
proximation is not in all respects valid. Using the improved potential evaluated at 
a Higgs boson mass of 35 GeV, one finds that the nucleation occurs at a temper- 
ature given by e ~ I where e = {T'^ — T'^)/{T^ ~T^)- Because the expansion of the 
Universe is so slow at this temperature, a typical bubble grows to a macroscopic size 
before colliding with other bubbles. The fractional change in the temperature of 
the Universe during the period of expansion is found to be roughly 5T/T ~ 10^^. 
It is thus a good approximation to ignore the expansion of the Universe in our 
calculations. 

In the first scenarios proposed for the formation of the asymmetry, baryon 
number was produced in the bubble wall. ' This mechanism, at best, is not 
terribly efficient, because the baryon number violating processes turn off rapidly as 
the scalar field expectation value turns on. The resulting asymmetry is sensitive to 
the speed and thickness of the bubble. The most effective scenarios for electroweak 
baryogenesis have the baryons produced in front of the wall, in the symmetric 
phase. ^^'^^ In this picture, scattering, for example, of top quarks from the bubble 



* We choose this value here and throughout this section because it represents roughly the 
largest value that would be compatible with baryogenesis. We emphasize that the minimal 
standard model is being taken here as a toy model. The strength of CP violation is of course 
another issue of importance for baryogenesis, but of no importance for the propagation of 
the bubble wall. 
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wall leads to an asymmetry in left vs. right-handed top quarks in a region near 
the wall. This asymmetry, resulting from an asymmetry between reflection and 
transmission of different quark helicities at the wall, biases the rate of baryon 
number violation in the region in front of the wall; the resulting value of n;,/n^ can 
be as large as 10~^. However, the authors of Ref. 15 assumed that the wall was 
rather thin, with a thickness of order T~^. For thicker walls, the asymmetry goes 
rapidly to zero. This can easily be understood. In order to have an asymmetry 
in reflection coefficients, the top quarks must have enough energy to pass through 
the wall. For mt ~ 120 GeV, this means typically the energy must be greater than 
about T/2. If the wall is very thick compared to this scale, the motion of the top 
quarks is to a good approximation semiclassical, and the reflection coefficient is 
exponentially suppressed. The analyses of other authors also exhibit sensitivity to 
the wall shape and velocity. 

Clearly, then, it is important to understand how the bubble propagates 
after its initial formation. A complete description of the wall evolution is rather 
complicated. We will see, however, that in certain limits it is not too difficult to 
determine how the velocity and thickness of the wall depend on the underlying 
model parameters. 

We begin with a history of calculations of the bubble velocity. In the 1970's 
Coleman^^ studied the zero temperature case and found that the bubbles quickly 
accelerated to the speed of light. At high temperatures however the Universe is 
filled with a plasma and one expects that it is important to consider the effects 
of this plasma on the propagation of the wall. There have been a number of 
studies of the kinematics of the process. "^^'^^ We will assume in this analysis that 
a steady state is reached after a sufficient time such that the wall propagates with 
a constant terminal velocity. In this case we can boost to the wall frame, in which 
the plasma on either side of the wall is characterized by a mean velocity 'yin(out) s-'^d 
temperature Tijj(out)- Applying conservation of energy-momentum to the plasma 
far on either side of the wall, one obtains two equations relating the parameters 
above (from T^^ and T*^^). The result of this analysis is that essentially two 
different types of propagation are possible: detonations, characterized by vin < f out 
and Tin > ^out; and deflagrations characterized by fin > fout and Tjn < Tout- 

In order to fully determine the velocities and temperatures, one needs ad- 
ditional dynamical information. One attains this by studying the interaction of 
the plasma with the wall. In the present case, the interaction is provided by the 
increase in mass of the various particles as they traverse the wall. It is clear that 
the most important particles here are the heaviest, namely the top quark and the 



f It may be possible to construct extended models for which this statement is relaxed. '^'^ 
-k The label 'in' refers to the broken phase, and 'out' to the symmetric phase. 
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W and Z . The lighter degrees of freedom essentially pass through the wall un- 
touched; however these degrees of freedom constitute approximately 80% of the 
plasma and are important in thermalizing the distributions of heavy particles. Be- 
cause of the preponderance of light species, it is natural to suppose that the changes 
in temperature and velocity across the wall are small. This assumption allows us 
to write a simple equation for the mean velocity of the wall. The smallness of the 
variations in velocity and temperature will be verified a posteriori. With these 
assumptions, the wall velocity is computed by summing up all of the forces on the 
wall and then solving for the velocity. In our discussion, we will assume that the 
velocity is non-relativistic and expand the force in small velocity. By looking at 
energy-momentum conservation in the wall frame, we can write 

d,T^' = J2^^'^ (2.1) 



where the left hand side contains the (zero temperature) energy-momentum tensor 
of the scalar field, and the right hand side is the sum of forces from the i species. 
At zero velocity, as we will see below, the right hand side is just the difference in 
pressure of the plasmas on either side of the wall, and thus the velocity independent 
part of the force is just the difference in the effective potential on either side of the 
wall. 

There have been a number of attempts in the past to compute the wall 
velocity. In Ref. 34 a simple formula for the wall velocity was given, based on a 
semiclassical picture in which a species of particle gains a large mass M ^ T as 
it passes through the wall. Balancing the force on the wall due to these particles 
with the pressure difference between the two phases gives a relation of the form 

V = ^^ (2^2) 

Ap 

where p is the pressure and p is the internal energy. 

More recently, Turok has argued"^^ that this type of analysis is incorrect. 
He suggests that reflection of particles from the wall does not slow the wall at all. 
The crucial step in Turok's original analysis is the assumption that equilibrium 
distributions appropriate to the value of the mass are maintained locally. Under 



I This assumption is not operationally crucial, but is merely convenient. 
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this assumption, the force on the wall can be written 



dz I „ n[k,x) V 



z 
z' 



(27r)3 ^ ' ^ " A^ 

, dm^ f Sk 1 -, (2.3) 

dz J {2nf2ko ^ ' 

= velocity — independent. 

The distribution n is taken to be the boosted Bose-Einstein (or Fermi-Dirac) dis- 
tribution. The second line follows by writing the change in momentum in terms of 
the change in mass as the particle traverses the wall (using energy conservation). 
The last line follows since the integral over momentum involves Lorentz-invariant 
quantities only, and thus the velocity dependence may be removed by a change of 
variables. With this analysis one concludes that the wall accelerates relentlessly, 
at least until some other effect takes over. This analysis is indeed correct; it is the 
assumption of equilibrium that is physically unmotivated. We will see shortly that 
the most important effects arise because of small, velocity-dependent departures 
from equilibrium. 

Before going on to the more realistic case hinted at in the last paragraph, we 
will improve on the analysis of Ref. 34; in the course of doing this we will identify 
some of the shortcomings of this analysis. 

2. 1. Thin Wall Analysis 

We will now make the following assumptions: we suppose that the individual 
particles of the plasma, while interacting with the wall, do not interact among 
themselves. We thus may follow individual particles as they traverse the wall and 
compute the force accordingly. This assumption obviously corresponds to assuming 
that any relevant mean free path is long compared with a scale set by the size of 
the wall. This, it will turn out, is an implicit assumption in the estimate of Ref. 34. 
However, we can improve on that analysis by noting that in reality there are no 
particles of relevance with mass much larger than the temperature. Depending on 
the momentum, a certain fraction will be reflected from the wall, but a significant 
fraction is also able to pass through the wall in either direction. If we then make the 
simple assumption that the incoming particles have a roughly thermal distribution, 
we can compute the total force on the wall by carefully adding up the contributions 
from each type of event noted above. The details of this calculation are given in 
Ref. 2, and we will simply quote the result here. If we write the force as 



Q = AVt + v ^£i + 0{ 



v') 
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we obtain, in the thin wall case, 



Si = p{o,T) - p{mi,T) - 



471-2 



E no{E)dE. 



(2.4) 



We note that in the limit m -^ oo, the result of Ref. 34 is recovered. If we expand 
this expression in powers of m/T, we find the following result. For top quarks, we 
have 



cthin ^ "-"/^ _,4 

and for W and Z bosons, we have 



In 






^7" 

2 


+ 


\t) 



[2.5) 



cthin 



30^ 



^ -ET<P' +..24 



.2 m2 



2,ni. ( lu '""'^ 



,2 j,2 



7V 4/1 "^10 ■ 



:2.6) 

If we ignore the top quark, the velocity turns out to be a pure number times a 
simple function of e. 



vr e 



61 



0.1 for m// = 35 GeV) 



valid for small e. Putting back in the top quark, we find that the velocity, in this 
approximation, is diminished. For mn = 35 GeV and mt = 120 GeV, one finds 
V ~ 0.07. 

We may now go back and check the validity of our assumptions. First, the 
working assumption of non-relativistic velocities clearly holds. Secondly, if we go 
back to the kinematics and compute the variations in velocity and temperature 
consistent with energy-momentum conservation we find 



6v 



V 



+ 10" 



^--10- 

T 

This result is clearly consistent with our assumption of slowly varying temperature 
and velocity. Furthermore, the signs of these quantities are such that we identify 
the process as a (weak) deflagration. 



* For the sake of this discussion, we are assuming that the expansion is vahd for top quarks. 
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2. 2. Parameters of Interest 

In order to proceed, we now investigate various distance scales relevant at 
this phase transition. To get an idea of the thickness of the wall, we note that at 
Tc there is an exact kink solution known of the form = 0c/2[l + tanh(2z/5)] and 
5 characterizes the wall thickness. In our case, we find 

E 

The value quoted is for m^j = 35 GeV. From numerical calculations, this width 
is accurate at the actual nucleation temperature as well, and thus we see that the 
wall is rather thick. 

It is useful to compare this number with the mean free paths for various 
processes. In considering the properties of the bubble wall, the relevant mean 
free paths are those for particles which interact with the wall, i.e., principally top 
quarks, W^s and Z's. The processes with the shortest mean free paths are elastic 
scatterings. These exhibit the characteristic singularities of Coulomb scattering at 
small angles. What actually interests us, however, is the momentum and energy 
transfer in these collisions. This is a problem which has been extensively studied, 
and we can borrow the relevant results. One finds that the momentum loss per 
unit length due to scattering is in all cases much larger than that due to the wall. 
This result may be understood in an alternative way. The elastic scattering cross 
section diverges at small angles in empty space. In the plasma, we expect that this 
divergence is cut off, essentially by a temperature factor. Examining the expression 
for the elastic scattering cross section, one obtains an estimate for the mean free 
paths of order £ ~ AT~^ for quarks, and £ ~ 12 T~^ for Ws and Z's. We will use 
these estimates in what follows but it should be kept in mind that a more complete 
analysis of the mean free paths is desirable. 

From the above discussion we see that the mean free paths for thermalization 
processes are typically a significant fraction of the wall thickness. We thus expect 
that a more realistic analysis of the bubble wall velocity lies somewhere in between 
the thin-wall analysis given in the previous section and the equilibrium analysis of 
Ref. 35. There have been two attempts at this in the literature. The first^ builds 
on the thin-wall analysis, whereas the second"^^ perturbs around the equilibrium 



f We make the assumption that the thickness of the waU remains constant throughout the 
expansion. 

X Also of interest are mean free paths for processes that change some approximately conserved 
quantum number. It has been speculated^ that a "snowplow" effect may occur in front of the 
wall which may further act to decrease the wall velocity. Estimates given elsewhere indicate 
that this effect is at most of approximately the same strength as the effects discussed herein. 
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picture. Both approaches show quahtative agreement that the wall velocity is 
non-relativistic. We will briefly discuss these two approaches in the following two 
sections. 

2. 3. DLHLL Thick Wall Analysis 

To get an idea of how finite elastic scattering lengths affect the velocity of 
the bubble wall, we assume that particles propagate freely over distances of order 
a mean free path, L As shown in the previous section this is typically shorter than 
the thickness of the wall and so we view the bubble wall as a succession of slices 
with thickness of order £, and for each of these we repeat the thin wall analysis. 
Thus we make the simple assumption that the distributions are given roughly by 
equilibrium distributions appropriate to a given point within the wall, and then 
follow that distribution over a length of order L On general grounds, one expects 
that the result may be written in the following form: 

c c cthin I c cthin 
^ = ^hH + '-'/^/ ) 

where S are suppression factors dependent on L One expects them to have several 
limiting behaviours. When (1/5 -^ they should approach zero, reproducing the 
equilibrium analysis in which the velocity dependence goes away. When 5 < l, the 
thin wall analysis is recovered. Numerically, one finds that the suppression factors 
are not very sensitive to mt and m//, and are well fit by 



Using the values quoted above for the mean free paths, the equations for £ from 
the previous section and assuming rriH ~ 35 GeV {5 ~ 40 T~^) and mt ~ 120 
GeV, we find a velocity of about v ~ 0.2 . Fig. 6 illustrates the velocity for a range 
of top masses, for two values of the Higgs boson mass. Also included on this plot 
for comparison is the thin wall result for these Higgs boson masses. 

2. 4. LMT Thick Wall Analysis 

An alternative analysis of the thick wall velocity has recently been per- 
formed in Ref. 36. These authors study linear perturbations of the distribution 
function n, determined by the relativistic Boltzmann equation 

d 

PzdzU + mjFz-^n = C{n). (2.7) 

Here the force term is given by the change in mass of the particles in the wall 

Fz = -l/2kodzm'^. 
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Let us now make the relaxation time approximation, whereby the colhsion 
term is replaced by C = 5n/ 1 where d, is again the relevant mean free path. One 
can then solve for 5n = n — no in terms of the equilibrium distribution function 
and the velocity of the wall. Balancing the forces on the wall one obtains the result 
of Turok, with an additional velocity dependent term 

/SVt = IV- — T(f)^. 
n 

Solving for the velocity, one obtains 

71 S e 

'-fV ^ c 



6 £ 1 + e' 

where c is of order unity, a result which obviously agrees very well with our result 
for very small i/6. 

The authors of Ref. 36 however perform a more detailed analysis, whereby 
they solve Eq. (2.7) directly without making the relaxation time approximation. 
The result of this analysis differs from the above in some respects, but within the 
range of parameters of interest agrees qualitatively with the above result. One 
must conclude from these analyses that the simple analysis of Ref. 2 is equivalent 
to the relaxation time approximation, but extends the analysis beyond very small 
i/6. The conclusion in all cases is that the velocity of the bubble walls is non- 
relativistic. 

3. Conclusions and Acknowledgments 

In this talk, we have discussed the recently improved understanding of the 
effective potential describing the electroweak phase transition of the minimal stan- 
dard model, and outlined several ways of understanding the resummation necessary 
to improve perturbation theory. We have also discussed recent calculations of bub- 
ble wall properties, understood how they relate to one another, and noted that all 
calculations are in qualitative agreement that the wall velocity is non-relativistic. 

I would like to thank Patrick Huet, Michael Dine and Andrei Linde for their 
collaboration on these projects and for many fruitful and interesting discussions. 
In addition, Andrei Linde was helpful in clarifying the historical record. I also wish 
to thank L. McLerran for discussion of his work. 
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FIGURE CAPTIONS 

1) Infrared divergent ring diagrams. 

2) Schwinger-Dyson equation. 

3) Tadpole graphs through two-loop order in Schwinger-Dyson approach. 

4) One-loop diagram in infrared effective theory, showing expansion in terms of 
ring diagrams. 

5) Typical electroweak potential showing hypothetical region of validity of ex- 
pansion. 

6) Plot of velocity versus mt for several values of mjj for thick and thin wall 
approximations. 
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